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Abstract: We discuss the modified gravity which includes negative and positive 
powers of the curvature and which provides the gravitational dark energy. It is shown 
that in GR plus the term containing negative power of the curvature the cosmic 
speed-up may be achieved, while the effective phantom phase (with w less than — 1) 
follows when such term contains the fractional positive power of the curvature. The 
minimal coupling with matter makes the situation more interesting: even 1/R theory 
coupled with the usual ideal fluid may describe the (effective phantom) dark energy. 
The account of R 2 term (consistent modified gravity) may help to escape of cosmic 
doomsday. 



*Also at TSPU, Tomsk, Russia and at IFT, UNESP, Sao Paulo (temporary) 



Contents 

|l|. Introduction [l] 

|2|. Sudden future singularity in modified gravity [2] 

|3|. Modified gravity coupled with matter [6| 

|3]. Discussion 



1. Introduction 

The interpretation of the very recent observational data indicates that current uni- 
verse is flat and is in accelerating expansion which started about five billion years 
ago. The important characteristic of the (accelerated) expansion is so-called "equa- 
tion of state" parameter w (for dominant energy density contribution) which is the 
relation between the pressure and the energy density of the universe. It is estimated 
that about 70 percent of universe energy density is composed of some mysterious 
effective fluid (dark energy which rules the universe) with negative pressure and w 
being close to —1. Despite the number of efforts we are still far away from the 
theoretical understanding of dark energy and its origin. 

The very promising approach to dark energy is related with the (phenomenolog- 
ical) modifications of Einstein gravity in such a way, that it would give the gravita- 
tional alternative to dark energy. Conceptually, it looks very attractive as then the 
presence of dark energy is only the consequence of the universe expansion. One such 
model containing 1/R term (which may originate from M-theory |2[]) was proposed 
in ref.|I| as gravitational alternative for dark energy. Modified gravity describes the 
accelerated expansion but it contains number of instabilities @. Nevertheless, fur- 
ther modification of the model by i? 2 -term [|J] or lni?-term || (see also ||) leads 
to consistent modified gravity which may pass solar system tests and is free of the 
instabilities. 1 One may consider other generalizations like ones including the positive 
(negative) powers of Ricci tensor squared ||, coupling of f(R) theory with scalar 
[§> 0) multidimensional 1/R theory [If | or even more extravagant (non-symmetric) 



1 It is interesting that Palatini formulation of modified gravity leads to physically different theory 
(for very recent discussion and list of refs., see 0) which seems to be free of (some) instabilities 
too. 
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gravity (for recent discussion, see |I2| ). Definitely, various predictions of consistent 
modified gravity should be tested. In its own turn, these tests may suggest further 
modifications giving true description of an observable universe. 

The present work is devoted to further study of the properties of modified grav- 
ity which contains positive and negative powers of the curvature. It is demonstrated 
that (effective phantom) dark energy cosmological solutions for the model naturally 
contain the finite-time, sudden singularity in the future (section 2). However, the 
consistent modified gravity with negative (or positive, fractional) powers of the cur- 
vature as well as with R 2 term has more stable future history. The cosmic doomsday 
does not occur there, rather the universe ends up in deSitter phase. In the third 
section it is considered the coupling of modified gravity with the usual matter. The 
very interesting observation shows that in such unified framework it is easier to re- 
alize various types of effectively dark energy universe. In summary, some outlook is 
given. 



2. Sudden future singularity in modified gravity 

Let us start from the general model of gravity depending only on curvature: 

S = -J d 4 x^f(R) . (2.1) 



Here f(R) can be an arbitrary function. It has been known [13] that such modified 
gravity theories can be rewritten into scalar-tensor form via the conformal trans- 
formation. Depending on the form of the function f the scalar-tensor theory may 
contain ghost-like term (with negative kinetic energy). It is quite remarkable that 
f(R) gravity of special form as we show below from this big class of theories may 
find the interesting applications as candidate to describe dark energy universe and 
its acceleration. In a sense, it is return of somehow forgotten generalized gravity 
motivated by recent astrophysical data. 

By introducing the auxilliary fields A, B, one may rewrite the action ( |2.1| ) as 
S = \ J d A Xy/—g {B (R — A) + f(A)}. Using equation of motion and expressing 
B in terms of A, one arrives at the Jordan frame action. Using the conformal 
transformation g^ u — > e^g^ with a = — In f'(A) (see also[TJ]]), we obtain the Einstein 
frame action [|]]: 

Se = -^J <?xV=g (« - \<ra r od^ - V(a)) , (2.2) 

V(a) = e° g (e~) - ef/ (g (e~) ) = - . (2.3) 

Note that two such theories in these frames are mathematically equivalent. However, 
physics seems to be different. For instance, in Einstein frame the matter does not 
freely fall along the geodesies which is well-established fact. 
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As an interesting example the action of (large distances) modified gravity may 
be taken in the following form 

S = -L J d A x^ (R - 1 RT n ) . (2.4) 

Here 7 is (an extremely small) coupling constant and n is some number. Since the 
function f(A) and the scalar field a are 

f{A) = A--fA- n , a = -\n{l + ni A- n - 1 ) (2.5) 

the potential is given by 

n 



l + wyA-™- 1 ) V n 7 

When curvature (~ A) is small and n > — 1 and n 7^ 0, the potential behaves as 
an exponential function V ~ ^ (l + i) ^4™+ 2 ^ (1 _)_ I) (77,7)"+! e"+ T<J . On the other 

1 1 n 

side, when the curvature is large, it follows V ~ r y(n + l)A~ n ~ (nj) n+l (— er) "+ 1 . 

Since V^'(A) = _ 7W (" +1 ) A iLJ^j^Z^ i ; \/(A) has only one extremum and the 

1 

extremum is given at A = {(n + 2)7} n+1 . 

The FRW universe metric in the Einstein frame is chosen as ds 2 E = —dt\ + 
a s(*s) Xli=i (<^ 1 ) • W the curvature is small, the solution of equation of motion 

3(n+l) 2 

IS CLE ~ ^ n+2) " . The FRW universe metric in Jordan frame is ds 2 = —dt 2 + 
a 2 (t) 5^i =1 (<ir J ) , where the variables in the Einstein frame and in the physical Jor- 
dan frame are related with each other by t = J e^dtE, a = ei cle- This gives t ~ t E 
and 

(n+i)(2n+i> Qn 2 + 7n — 1 . „. 

a ~ t ™+ 2 , w = — . (2.7) 

3(n + l)(2n + l) 



1 

n+2 



The first important consequence of above Eq. (|2.7|) is that there is possibility of 
accelerated expansion for some choices of n (effective quintessence). In fact if n > 
~ 1 "^ v ^ or — 1 < n < — |, we find w < — | and > 0. 2 In other words, modified 
gravity presents the gravitational alternative for dark energy with the possibility of 
cosmic speed-up. The corresponding analysis was given in detail in ref.pl. 

When — 1 < n < — ~, it follows w < — 1. If w < — 1, the universe is shrinking 
in the expression of a Q2.7| ). If we replace the direction of time by changing t by 



2 In the solution (2.7), the Hubble parameter H = has the form oi H = ^f-, where /iq is 
a constant of the unity order. There is an ambiguity how to choose t in H. One natural choice 
is to take t to be of the order of the age T of the present universe. Since T ~ 1.37 x 10 10 years 
~ (l0~ 33 eV) , we find H ~ 10 -33 eV, which is consistent with the observed value of H in the 
present universe: -^observed ~ 70kms _1 Mpc _1 ~ 10~ 33 cV. Then even if the original Lagrangian 
theory does not contain small parameter of the order of the Hubble parameter, such a small scale 
is naturally induced from the age of the universe in the power law expansion (2.7). 
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—t, the universe is expanding but t should be considered to be negative so that the 
scale factor a should be real. Then there appears a singularity at t — 0, where the 

2 

scale factor a diverges as a ~ (— t) 3 (™+ 1 ) . One may shift the origin of the time by 
further changing —t with t s — t. Hence, in the present universe, t should be less than 

2 

t a and it looks to appear the singularity at t = t s : a ~ (t s — t) 3 ( m+1 ). The future 
singularity may be called sudden or Big Rip singularity. We should note, however, 
the expression of w in fl2.7| ) could be correct only when curvature is small, as in the 
present universe. When t goes to t s , the curvatures become large as (t s —t)~ 2 and 
Eq. fl2.7Q becomes invalid. As the curvature becomes large, the first Einstein-Hilbert 
term in (|2.4j ) dominates. Then as will be shown later, the singularity is moderated 
and if we include usual matter, the singularity will not appear eventually. Since 
w + 1 = ^J+rjl^n+i) > ft f°U° ws w < —1 when — 1 < n < — \ (the effective phantom 
phase for fractional positive power of the curvature). Here it is assumed n > —1, so 
that the Einstein term dominates in ( |2.5| ) when the curvature is small. We should 
also note when n > 0, as ^J+i){2n+i) > 0, w is always greater than —1. Then the 
theory with negative power of the curvature like -^-gravity does not produce effective 
phantom although it may produce effective quintessence. 

The results ( |2.7| ) are valid when the curvature is small but near the Big Rip sin- 
gularity, the curvature becomes large and (|2.7|) is not valid. The qualitative behavior 
when the curvature is large can be found from the potential. 

The qualitative behavior of the potential when — 1 < n < — | and 7 > is given 
in Figure |l|a). In order that a is real, however, Eq. fl2.5Q tells R ~ A > (—wy)™^ 1 . 
Then the curvature cannot be small and the expressions fl2.7Q are not valid. When 
A < (— n7)"+ T , instead of ( |2.5| ), one may define a = — In (— 1 — n / -fA~ n ~ 1 ). As 
shown in ||, however, the anti-gravity appears in this case and instead of ( |2.2Q , 
we obtain Se — ^ / d 4 Xy/^g + \g po 'd p ad a a — V{o)\ with potential V{a) = 
—e a g(e~ a ) — e 2a f (g (e™ ")). Hence, when — 1 < n < — \ and 7 > 0, the region 
A < (— 727) "+ 1 is not physical. Then we should assume 7 < 0. In case —1 < n < 
and 7 < 0, there does not appear the extremum in the potential if curvature is 
positive. The qualitative behavior of the potential is given in Figure |I]b). The 
Eq.( |2.6|) shows that the potential is negative and is unbounded below when A is 
large. 

In order to consider the region where the curvature is large for the case — 1 < 
n < — I and 7 < 0, the following potential is taken 

V = -V (-a) a , V = —( ni )— >0, a =JL-<-l. (2.8) 

n n + l 

The cr-equation of motion and the FRW equation in the Einstein frame are = 
~ 3 (If + 3H lt) - V '( a ) and % H e = I (ft) + V ( a )- A consistent solution is 
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Ci C 2 V(A) 

a) b) 
Figure 1: a) A qualitative behavior of the potential with — 1 < n < — | and 7 > 0. Here 

1 1 1 

C\ = (—nj) "+ 1 and C2 = (— (n + 2) 7) "+ 1 In the region with A < (—ivy) n+1 , there appears 

the anti-gravity phase, b). A qualitative behavior of the potential with — 1 < n < — ^ and 

7 < 0. The potential is monotonically decreasing. 



given by 



1 



H E = 0, „ = ( t , E - t E f ,^ s (™>Y- llis * = 2 ±±±> . (2.9) 

\op / 2 — a n + 2 

Here £ s ^ is a constant of integration. Since — 1 < n < — ^, it follows < /? < | 
. The spacetime in the Einstein frame is flat and the scale factor a E = ao is a 
constant which is not the case in the physical Jordan frame, where they are finite 
even at t E = t sE although there is a cut singularity there. The cut singularity, 
however, makes the physical curvature divergent. In fact, when t E ~ t sE , t is given 
by t = t s + t E — t S E with a constant of the integration t s . In the Jordan frame, the 
scale factor a is given by a{t) ~ aoe~^ ( -* s_ ^' 3 and the Hubble parameter H = 
behaves as H = |f = = ^ (t fljB - " x e"f ~ ^(t.-t)^ 1 . Since 

— 1 < /3 — 1 < — 3 < 0, diverges at t = t s as well as the scalar curvature 
R = Q^- + 12H 2 . In case of the Big Rip singularity, H behaves as H ~ ^3^. Hence, 
the behavior of H is moderated. In case of the Barrow model |l5j (see also []16, p~Tf] ) , 
H behaves as H ~ h + /ij (t s — t) a with constant /iq and hi. As < a < 1, the 
behavior of H here is more singular than that in [15 . The singularity in the present 
case is moderated and a does not diverge. With the usual matter inclusion such 
a singularity will not appear. It is interesting to note that for the model [[OJ, the 
singularity is moderated (or even disappears) if the quantum corrections are taken 
into account |ll| . Note also that usual, infinite-time singularity is still possible even 
in consistent modified gravity (for earlier, related discussion of infinite time, future 
singularity in f(R) gravity, see jl9f). 

In accord with the proposal of ref . J| one may add a term proportional to R 2 to 
the action (|2~4l ) 

S = -^ J d 4 x^ (R - 7 iT n + r]R 2 ) . (2.10) 

In this situation the unified theory permits early time inflation as well as late time 
acceleration and does not contain the number of instabilities. The function f(A) (|2.5|) 
is modified to be f(A) = A — jA~ n + i]A 2 . We now assume — 1 < n < 7 < 0, and 
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r] > 0. Then a = - In (1 + nyA' 71 " 1 + 2r)A). Although l + nyA^ 1 + 2r)A > 0, since 

s = "7W = there is a branch point ' where s = or = 

0, at A = A = {KE+ih}^ or a = ao = -In (l + (n + 2) (^_)^ (n 7 )^) . The 
potential V(.A) has the following form: 

(1 +n 7 A~™- 1 + 2^A) 2 



The behavior when A is small is not changed from the case in ( |2.4| ). On the other 
hand, when A is large, V(A) goes to a constant: V(yl) — > y. The potential V(A) 

vanishes at A = Ai = |^±lh|~. Note ^ Q < ^ since < m < _ 7 . y'(A) 

has an extremum at A = A since V'(A) = i n( " +1)7 ^ 1+n7 ^-n-f+2r ? l"3 +2hA -• 

Since there is a branch point at a = cr ,if we start from the small curvature, the 
growth of the curvature stops at R — A , where a = Oq. In fact, at the branch 

d 2 V 



point, where f"(A) = 0, the mass m a of a becomes infinite since m a oc 

/"(A)<L4 djT) - f"(A) + /'(A) + /'(A)2 ^ +°°- iNote alS0 J l U J < U wnen 

A < Aq. Then the growth of a is finished at a = a . Hence, adding R 2 term, 
there does not occur cosmic doomsday but the universe ends up in deSitter phase. 

The scale factor a is given by a(t) ~ aoe V 12 with a constant a . Note that the 
quantities in the Einstein frame are different from those in the Jordan frame by 
almost constant factor as ag ~ e~ or dt E ~ e~dt. This supports our point of view 
that cosmic doomsday in such theory does not occur because after the late time 
acceleration the universe starts new inflationary era. It is also important to stress 
that as the mass of a becomes very large, there is no problem about the equivalence 
principle since a cannot mediate the force. In other words, unlike to BD theory, such 
a model may pass the solar system test provided by VLBI experiment. 



3. Modified gravity coupled with matter 

It is very interesting that modified gravity which can be made consistent one [|J 
may help in the resolution of dark energy problem in various ways as it suggests 
gravitational alternative for dark energy. In particular, as we will show below it 
may give phantom dark energy without necessity to introduce the (negative kinetic 
energy) phantom scalar theory. In fact, the matter is taken to be the usual ideal 
fluid. 

We now consider the system of the modified gravity coupled with matter: 

S= [ d 4 x^{f(R)+L rn } . (3.1) 
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Here f(R) is an adequate function of the scalar curvature and L m is a matter La- 
grangian. Then the equation of the motion is given by = ^g ia> f(R) — R^f^R) — 
V^Vj, /'(.R) — g^yV 2 f'(R) + \T^ V . Again, the FRW spacetime is considered. The 
ideal fluid is taken as the matter with the constant w: p = wp. Then from the 
energy conservation law it follows p = poa~ 3 ^ 1+w ^ . In a some limit, strong cuvature 
or weak one, f(R) may behave as f(R) ~ foR a , with constant fo and a. An exact 
solution of the equation of motion is found to be 

I, _ 2a 
a = a t , h 



a 



3(1 + w) 



(-6/i + 12/^o) Q {(1 - 2a) (1 - a) - (2 - a)h } 



.(3.2) 



Po 

When a = 1, the result h = 37^3 in the Einstein gravity is reproduced. Note that 
stability issue should be carefully investigated here. However, even if the solution 
is instable the decay time could be very big due to the fact that coupling constant 
of modified gravity term is very small. From another side, when finite-time future 
singularity occurs it may be resolved by the account of quantum effects [|18|, || . 



The effective w e s may be defined by ho = 3 ( 1+w ff ) • By using fl3.2|) , one finds 

w efi = -H . (3.3) 

a 

Hence, if w is greater than —1 (effective quintessence or even usual ideal fluid with 
positive w), when a is negative, we obtain the effective phantom phase where w e ff is 
less than —1. This is different from the case of pure modified gravity. Moreover, when 
a > w + 1 (it can be even positive), Wes could be negative (for negative w). Hence, it 
follows that modified gravity minimally coupled with usual (or quintessence) matter 
may reproduce quintessence (or phantom) evolution phase for dark energy universe 
in an easier way than without such coupling. 

One may now take f(R) as in ( |2.1(J| ). When the cuvature is small, the second 
term becomes dominant and one may identify fo = — \ and a = —n. Then from 
(P-3|) , it follows w e ff = —1 — _ Hence, if n > 0, we have an effective phantom 
even if w > — 1. Usually the phantom generates the Big Rip singularity. However, 
near the Big Rip singularity, the curvature becomes large and the last term becomes 
dominant. In this case a = 2 and w e g = — . Then if w > —1, it follows 
w e g > — 1, which prevents the Big Rip singularity (makes phantom phase transient) 
as is described in the previous section. To conclude, it looks quite promising that 
modest modification of General Relativity coupled to ideal fluid matter leads to 
effective dark energy universe in the very natural way. 

4. Discussion 

In summary, the gravity theory with negative (like 1/R) and positive (quadratic) 
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powers of scalar curvature shows the number of features which are desireable to 
explain the accelerating dark energy universe: 

1. It passes solar system tests (VLBI experiment). Thanks to higher derivative R 2 
term, the gravitationally bound objects (like Sun or Earth) are stable as well as in 
GR. Of course, some fine-tuning of coupling constants is necessary as is shown in 
||]. (It was mentioned already in ref.[]I] that the coefficient of l/i?-term should be 
extremely small.) 

2. Newtonian limit is recovered just at the above values for coupling constants. 

3. The gravitational dark energy dominance is explained simply by the universe 
expansion. Moreover, when modified gravity is coupled with usual matter it is easier 
(less deviations from GR) to get the effective (phantom or quintessence) dark energy 
universe regime. 

4. The presence of R 2 term in consistent modified gravity prevents the development 
of the cosmic doomsday. 

Thus, the consistent modified gravity remains to be viable candidate for the 
explanation of dark energy as the gravitational phenomenon. Nevertheless, only 
future, more precise astrophysical/gravitational data will prove if it is the time for 
new gravitational physics to enter the game. 
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